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摘 要 : 作物 水 分 生产 率 的 准确 和 定量 化 评价 是 提高 干旱 区 作物 产量 的 基础 ,对 缓解 水 资源 短缺 和 
农业 可 持续 发 展 上 共有 重要 意义 。 以 塔里木 盆地 北岸 的 渭 干 河 - 库 车 河 绿 洲 ( 渭 库 绿洲 ) 为 典型 区 
域 , 基 于 Google Earth Engine(GEE) 云 平台 ,通过 建立 2009 一 2020 年 流域 SEBAL E RAWEA, 
棉花 分 布 识别 模型 及 估 产 模型 ,对 流域 棉花 水 分 生产 率 进行 评价 。 结 果 表 明 :(1) 渭 库 绿洲 棉花 产 
量 从 2009 年 的 1610.10 kg-hm?34 K 4! 2020 4 AY 1855.05 kg:hm”, 增 长 率 为 13.20%, 棉 花 种 植 面积 
永年 向 绿洲 边缘 延伸 ,棉花 产量 重心 整体 自 西向 东 移 动 2485 mo(2) 棉花 生长 期 2009 年 蒸 散 发 均 
值 为 686.80 mm,2020 年 为 738.66 mm, 整 体 呈 上 升 趋势 ,其 增长 率 为 7.02% ,棉花 生长 期 蒸 散 发 最 
大 值 为 花 铃 期 和 吐 你 期 , 蒸 散发 较 高 值 主要 分 布 在 绿洲 内 部 与 塔里木 河北 岸 讽 缘 。(3) 2009 年 水 分 
生产 率 均值 为 0.21 kg*m’,2020 年 均值 为 0.25 kg.m?,12 a 间 水 分 生产 率 均 值 增长 率 为 16%。 在 空 
间 上 , 涓 库 绿洲 水 分 生产 率 重 心 在 红旗 镇 自 东北 向 西南 移动 1832 m, 年 均 移动 速度 为 152.67 ma 
绿洲 棉花 水 分 生产 率 呈 现 东西 方向 大 于 南北 方向 扩张 趋势 ,空间 分 布 方向 趋势 增强 ,空间 格局 趋 
向 集聚 化 。(4) 12 a 间 产量 的 增长 速度 超过 了 蒸 散发 的 上 升 速度 ,促使 水 分 生产 率 提 高 。 其 次 ,水 
分 生产 率 与 棉花 种 植 面积 和 合理 的 水 量 灌溉 技术 密切 相关 ,水 分 生产 率 高 值 主要 分 布 于 沙 雅 县 新 
晨 农 场 和 新 和 县 桑 塔 木 农场 ,由 于 农场 规模 化 种 植 和 集约 化 管理 ,促进 了 棉花 增产 、 农 业 水 资源 的 
稳定 分 配 和 高 效 利用 。 

关 键 词 : Google Earth Engine(GEE); BRR; 棉花 提取 ; 棉花 产量 ; 
文章 编号 : 1000 - 6060(2023)10 - 1632 - 11(1632 ~ 1642) 


水 分 生产 率 ; 渭 库 绿洲 


栅 花 作为 最 重要 的 经 济 作物 ,在 我 国 新 疆 广泛 水 分 生产 率 是 指 单位 水 资源 量 在 一 定 的 作物 


种 植 ,2022 年 新 疆 棉花 种 植 面积 为 249.69x10’ hm’, 
总 产 达 539.1x10t, 棉 花 种 植 面积 .单产 、 总 产 连续 
28 a 位 居 全 国 首位 "1, 在 新 疆 经 济 发 展 中 占有 举 足 
轻重 的 地 位 。 然 而 长 期 以 来 干旱 少 十 和 缺 水 极 大 
的 限制 了 该 地 区 覃 花 健 康生 长 ,尤其 近年 来 ,区 域 
内 生活 、 工 业 和 生态 用 水 占 比 逐渐 增加 ,农业 用 水 
占 比 及 可 用 水 量 在 不 断 减 小 ” ,导致 新 疆 棉 花生 产 
面临 巨大 挑战 。 因 此 ,解决 水 资源 短缺 与 提高 作物 
产量 ,以 及 提高 农业 用 水 效率 ,已 成 为 新 疆 尤 其 是 
南 疆 农业 发 展 中 面临 的 关键 难题 所 在 。 
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品种 和 耕作 栽培 条 件 下 所 获得 的 产量 ,是 作物 灌溉 
管理 和 水 资源 优化 配置 的 主要 依据 和 核心 ”“。 定 量 
评价 区 域 作物 水 分 生产 率 是 缓解 干旱 区 水 资源 短 
缺 和 生态 环境 危机 的 重要 前 提 。 以 往 水 分 生产 率 
人 研究 基于 田间 实验 观测 和 采用 行政 单元 统计 ,不 能 
准确 表达 空间 分 异 和 区 域 特征 的 需求 。 随 着 3S 技 
术 和 Google Earth Engine(GEE ) 等 平台 的 发 展 ,遥感 
在 农业 领域 的 应 用 日 益 广泛 汪 , 到 感 估算 作物 耗 水 
量 理论 和 方法 日 益 完善 ,为 区 域 作物 识别 .产量 及 
耗 水 量 估算 等 提供 了 可 能 ,其 中 MODIS 遥感 数据 具 
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A es EY Ti] OP ES JE EE A i EE Se Fh 
用 MODIS 数 据 对 小 麦 、 玉 米 和 大 麦 等 水 分 生产 率 分 
析 了 时 空 变 化 ,并 探究 了 灌溉 系数 与 影响 因素 。 然 
而 ,MODIS 数据 的 空间 分 辨 率 有 限 ,对 于 破碎 化 耕 
地 难以 达到 精细 要 求 %。Landsat 影像 有 更 高 的 空 
间 分 辨 率 , 可 以 满足 长 时 间 序 列 的 耕地 作物 识别 精 
度 ,Zhou 等 "和 Ghorbanpour 等 号 利用 遥感 数据 佑 
算 了 冬小麦 和 夏 玉 米 等 水 分 生产 率 ,进一步 对 管理 
尺度 评价 ,但 以 上 人 研究 中 缺少 对 棉花 作物 水 分 生产 
率 的 分 析 , 虽 然 Mauget"”、Chen $E" Thorp 等 由 和 
Himanshu 等 “使 用 田间 控制 灌溉 的 方法 对 小 型 棉 
田 水 分 生产 率 估算 ,但 研究 集中 于 小 区 域 ,缺乏 较 
大 空间 尺度 长 时 间 序 列 的 研究 ,并 无 法 满足 棉花 
水 分 生产 率 时 空 变化 的 分 析 视 角 。 

鉴于 此 ,本 文 利用 GEE 平 台 的 优势 构建 棉花 识 
别 模型 . 估 产 模型 和 燕 散 发 模型 ,对 渭 干 河 - 库 车 河 
绿洲 (简称 渭 库 绿洲 ) 棉 花 水 分 生产 率 进行 探究 , 揭 
示 棉 花 种 植 和 产量 时 空 分 布 棉花 生育 期 内 耗 水 量 
的 时 空 变化 规律 ,在 此 基础 上 对 棉花 水 分 生产 率 进 
行 定量 评价 ,可 为 干旱 区 棉花 生产 提供 科学 管理 和 
技术 参考 。 


1 研究 区 概况 


渭 库 绿洲 (82°00’ ~ 83°60'E,41°00’ ~ 41°85'N) 
位 于 新 疆 南 部 ,地 处 塔里木 盆地 和 塔克拉玛干 沙漠 
北 缘 , 是 典型 而 完整 的 山 前 冲 洪 积 扇 平 原 ,隶属 阿 
克 苏 地 区 的 库 车 、 沙 雅 和 新 和 3 个 县 级 行政 区 (图 
1)。 光 热 资源 丰富 ,降水 稀少 ,蒸发 强烈 ,气温 差异 
显著 ,年 均 气 温 9.7~12.9 %C, 年 降水 量 28.9~194.7 mm, 
年 均 蒸发 量 2000.7~2092.0 mm" 1。 渭 干 河 成 为 流域 
内 生态 、 农 业 、 人 民生 产生 活用 水 的 主要 来 源 "。 
渭 库 绿洲 是 干旱 区 绿洲 农业 的 典型 代表 ,以 种 植 棉 
花 . 小 麦 作物 为 主 。 


2 数据 与 方法 


2.1 研究 数据 与 处 理 

为 保证 研究 区 该 年 有 充足 的 影像 估算 棉花 全 
育 期 的 蒸 散发 ,结合 2014 年 国家 政策 对 新 疆 棉花 产 
量 和 价格 的 调整 ,选用 前 后 年 份 的 数据 , 即 2009 一 
2020 年 Landsat 影 像 逐 年 6 一 7 月 数据 用 于 棉花 提取 
和 产量 估算 ;其 次 筛选 2009 年 ( 含 2008 年 5 月 8 日 、 


图 1 研究 区 概况 
Fig. 1 Overview of the study area 


2008 年 6 月 9 日 ) 和 2020 年 ( 含 2019 年 4 月 5 日 ) 棉 
花生 长 期 (4 一 10 月 ) 有 影像 14 景 ,结合 使 用 气候 再 分 
析 数 据 ERA5 、 数 字 高 程 数 据 SRTM EHTA 
感 反 演 ; 同 时 使 用 蒸 散发 数据 MOD16A2 和 气象 数 
据 验证 蒸 散发 的 精度 ,Landsat 数 据 基本 信息 见 表 1。 


表 1 研究 数据 
Tab.1 List ofthe data used in this study 


空间 分 ”波段 影像 

名 称 Wem 数量 成 像 时 间 ( 年 -月 -日 ) 数量 
Landsat 8 30 13 2013-06-23、2014-07-08 14 
2015=07=15、2016=07=17 
2017-07-04 .2018-07-23 
2019-07-26 .2019-04-05 
2020-05-25 .2020-06-26 
2020-07-28 .2020-08-29 
2020-09-14 ,2020-10-16 
2008-05-08 ,2008-06-09 9 
2009-04-25 2009-07-14 
2010-08-31 ,2009-09-31 
2009-10-18 ,2010-07-17 
2011-07-04 
2012-05-11 .2012-07-14 2 


Landsat_5S 30 13 


Landsat_7 30 11 


实地 调研 数据 来 自 研 究 区 历年 农作物 调查 数 
据 ,2009 年 至 今 每 年 一 次 ,获取 30 个 样 方 , 样 方 大 小 
为 500 mx500 m, 样 方 包括 作物 类 型 .面积 .长 势 .水 
分 温度 等 ; 另 历年 采集 各 类 作物 随机 样本 点 若干 。 

气象 数据 来 源 于 国家 气象 信息 中 心 (http://data. 
cma.cn) 发 布 的 渭 库 绿洲 内 3 个 气象 站 ( 沙 雅 站 、 库 
车 站 新 和 站 ) 的 逐日 平均 气温 、 最 低 、 最 高 气温 、 降 
雨 . 风 速 . 日 照 时 间 等 ;2009 一 2020 年 《阿克苏 统计 
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年 鉴 》 作 为 棉花 产量 估算 参 考 和 验证 数据 。 计算 公式 如 下 : 
遥感 数据 均 来 自 于 GEE 云 平 台 (https:Wearthen- y =y NDVI, (2) 
a NDVI 


gine.google.com ) ,利用 GEE 云 平 台 对 影像 数据 去 云 、 
裁剪 和 拼接 等 预 处 理 , 将 空间 分 辩 率 统一 重 采 样 为 
30 m, 然 后 使 用 该 平台 提供 的 JavaScript 语 言 和 消 数 
进行 数据 处 理 和 模型 构建 。 
2.2 研究 方法 

本 研究 根据 作物 水 分 生产 率 原 理 , 通 过 遥感 数 
据 等 数据 ,在 GEE 平 台 建 立 研 究 区 域 棉花 种 植 区 分 
布 识别 模型 . 估 产 模型 和 SFEBAL 遥感 蒸 散 发 模型 ， 
进而 构建 棉花 水 分 生产 率 估算 模型 。 同 时 结合 局 
部 空间 自 相关 标准 差 槛 圆 与 重心 迁移 分 布 等 地 理 
分 析 方 法 ,并 准确 和 定量 化 评价 研究 区 棉花 水 分 生 
产 率 时 空 分 布 特征 。 
2.2.1 棉花 水 分 生产 率 估算 水 分 生产 率 反 映 了 农 
业 水 量 投入 和 产 出 的 效率 ,是 衡量 农业 科学 和 合理 
用 水 的 综合 指标 , 即 作 物 经 济 产量 与 实际 蒸 散发 的 
比值 ,计算 公式 如 下 : 


22e 
veS T (1) 


式 中 : WP 为 作物 水 分 生产 率 (kg.m”);Y 为 作物 产 
Ht (kg hm) ,通过 建立 遥感 估 产 模型 估算 得 到 ; ET 
为 作物 的 蒸 散 发 (mm.d) ,通过 遥感 定量 反 演 获 
得 ,并 通过 气象 站 和 莹 散发 产品 数据 验证 。 
2.2.2 棉花 产量 估算 根据 预 处 理 后 的 Landsat 系列 
影像 ,选取 各 地 物 样本 点 的 纯净 像 元 ,利用 CEE 平 
台 构 建 主要 作物 生育 期 归 一 化 植被 指数 (NDVD"” 时 
间 序 列 曲 线 , 根 据 棉花 NDVI 特 征 , 获 取 棉 花 出 苗 期 
和 花 铃 期 的 影像 ,通过 时 间 特 征 ,结合 光谱 特征 AC 
理 特征 .地形 特征 等 构建 分 类 特征 。 将 提取 的 枫 花 
种 植 面积 与 统计 年 鉴 数据 进行 逐年 对 比 , 其 平均 相 
对 误差 均 小 于 5$%。 最 终 ,选取 分 类 精度 较 高 的 随 
机 森林 方法 进行 分 类 ,总体 精度 为 83% ,Kappa 系 数 
为 0.78。 因 此 ,本文 利用 随机 森林 分 类 方法 提取 研 
究 区 棉花 种 植 面积 。 

进一步 使 用 NDVI 指 数 进行 产量 估算 ,以 往 
研究 中 棉花 的 产量 与 NDVI 具 有 较 强 的 相关 性 ,并 
在 实际 得 到 充分 的 验证 了 3, 因此 本 研究 使 用 栅 格 
图 像 推算 ,将 县 级 产量 图 扩展 到 像 元 级 空间 分 布 ， 
结合 统计 年 鉴 中 的 棉花 经 济 产量 数据 ,得 到 研究 区 
3 个 县 级 区 划 的 棉花 产量 像 元 级 的 均值 ,然后 利用 
棉花 的 NDVI 指 数 扩展 到 单个 像 元 的 棉花 产量 ， 


式 中 : 7 为 单个 像素 级 的 棉花 产量 (kg hm ”); Yu 
为 县 级 棉花 平均 产量 ,从 统计 年 鉴 获得 ; NDVI, 为 
花 铃 期 NDVI 值 ; NDVI,,, 为 县 级 花 铃 期 平均 NDVI 
值 。 然 后 根据 数值 建立 棉花 产量 线性 关系 。 

2.2.3 MH AMR BRI SEBAL (Surface ener- 
gy balance algorithm for land ) 模 型 是 以 陆 面 能 量 平衡 
为 理论 基础 ,通过 遥感 影像 和 气象 参数 对 各 地 表 分 
量 进 行 估算 ,得 到 区 域 的 实际 蒸 散 发 。 地 表 能 量 平 
衡 方程 为 ”: 


AET=R,-G-H (3) 
IP: R, 为 地 表 的 净 辐 射 通 量 (W.m”); G 为 土壤 
热 通 量 (W.m3); 克 为 下 垫 面 到 大 气 的 感 热 通 量 
(W-m”’); AET 为 潜 热 通 量 (W.m”), 其 中 4 为 潜 热 
蒸发 系数 , ET 为 蒸 散 发 (Wm”) ,在 计算 水 分 生产 
率 时 ,将 ET 单位 转换 为 mmd'。 
因为 卫星 过 境 时 所 观测 的 是 地 面 瞬 时 数据 ， 
此 把 瞬时 薰 散 发 扩展 到 日 蒸 散发 , 蒸 散 比 在 白天 变 
化 较 小 ,可 以 假定 白天 蒸发 比 保持 不 变 , 根 据 瞬时 
蒸 散 比 和 日 有 效 辐射 能 量 计算 出 日 莱 散 发 ,利用 日 
蒸 散 比 插值 法 ,得 到 卫星 过 境 日 前 后 一 段 时 间 的 逐 
HRA” ,然后 根据 各 日 的 有 效 净 辐 射 量 ,计算 
出 逐日 蒸 散发 : 


ET, = = RT, (4) 


式 中 : ET, AS d RAHA; ET, 为 卫星 过 境 时 i 
FEE AS A ACE o 

2009—2020 4F [i] ,由 于 缺乏 中 间 年 份 棉花 生长 
期 蒸 散 发 估算 的 连续 遥感 影像 ,本文 根据 棉花 生长 
物候 历 ,选取 起 始 和 最 终 2 a 的 棉花 生长 期 影像 估 
算 蒸 散发 。 以 研究 区 新 和 ORE . 库 车 3 个 气象 站 点 
作为 样 点 ,使 用 基于 气象 站 的 P-M 方 法 与 SEBAL 模 
型 对 3 个 站 点 的 平均 蒸 散发 精度 验证 ,得 到 2009 
年 蔡 散 发 平均 绝对 误差 为 0.56 mm + d” 、2020 年 为 
0.27 mm d ,同时 利用 MODIS16A2 数据 8 d 产品 棉 
花生 长 期 累积 蒸 散发 精度 检验 ,与 SEBAL 模 型 反 演 
的 2009 年 蒸 散发 平均 绝对 误差 为 0.23 mm- d” 2020 
年 为 0.27 mm'.'d ,两 种 数据 都 具有 较 高 一 致 性 ,并 
取得 较为 满意 的 估算 结果 。 
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3 结果 与 分 析 
3.1 棉花 产量 时 空 分 析 

棉花 产量 的 时 空 变化 是 研究 水 分 生产 率 的 前 
提 , 将 遥感 估算 的 县 级 总 产量 与 统计 年 鉴 数 据 逐 年 
精度 对 比 ,其 平均 误差 小 于 8%。 如 图 2 所 示 ,2009 一 
2020 年 渭 库 绿洲 棉花 产量 呈 增 加 趋势 ,从 2009 年 
的 1610.10 ke-hm ”增长 到 2020 年 的 1855.05 kg-hm”, 
产量 增长 率 为 13.20%。 从 各 县 来 看 ,新 和 县 增长 
6.09% , 沙 雅 县 增长 12.729%5 , 库 车 市 增长 20.749% ,新 
和 县 的 产量 最 高 ,其 次 为 沙 雅 县 和 库 车 市 。 渭 库 
绿洲 产量 在 2010 年 最 低 (1520.04 kg+hm™) , 2019 4 
达到 最 高 (1984.80 kg + hm?) ,2014 年 出 现 低谷 
(1645.05 kg.hm?)。 由 于 2014 年 棉花 价格 调整 , 开 
始 大 面积 种 植 ,水 土 和 管理 措施 不 成 熟 等 原因 使 得 
棉花 产量 较 低 ,而 在 2015 年 之 后 趋 于 稳定 上 升 , 主 
要 由 于 棉花 种 植 面积 增加 导致 产量 逐年 上 升 汪 ,加 
之 政策 补贴 .水 肥 精 细 管 理 和 机 械 化 程度 持续 提 
高 ,对 棉花 产量 起 了 推动 作用 。 

采用 500 mx500 m 格 网 ,对 2009 一 2020 年 逐年 
棉花 产量 按 县 级 行政 进行 重心 迁移 规律 分 析 ,新 和 
县 产量 重心 在 玉 奇 喀 特 乡 和 渭 干 乡 之 间 移 动 (图 
3a)。 沙 雅 县 棉花 产量 重心 稳 居于 古 勒 巴 格 乡 ( 图 
3b) ,呈现 出 先 西 后 东 移动 趋势 。 库 车 市 产量 重心 
主要 位 于 墩 阔 坦 镇 、 阿 克 下 斯 塘 乡 、 齐 满 镇 (图 
3c)。 渭 库 绿 洲 产 量 重心 在 沙 雅 县 红旗 镇 和 努 尔 巴 
格 乡 之 间 波 动迁 移 ( 图 3d) ,12a 间 自 西向 东 移 动 
2485 m, 表 明 渭 库 绿 洲 棉花 产量 随 着 种 植 面积 增加 


何 旭 刚 等 : 基于 GEE 平 台 涓 库 绿洲 棉花 水 分 生产 率 遥 感 估算 


1635 


向 东西 方向 扩展 。 总 体 上 ,棉花 产量 各 时 间 点 的 标 
准 差 椭圆 长 轴 各 自 沿 绿洲 下 游 方 向 逐渐 移动 , 库 车 
市 和 沙 雅 县 椭圆 面积 大 于 新 和 县 ,新 和 县 总 体 趋势 
自 东 北向 西南 方向 扩展 , 沙 雅 县 自 西 向 东 趋 势 移 
动 , 库 车 市 则 自 西南 向 东北 方向 延伸 (图 3e)。 
3.2 棉花 生长 期 蒸 散 发 量 时 空 分 析 

棉花 生长 期 营 散 发 的 准确 估算 是 研究 水 分 生 
产 率 的 重要 环节 , 渭 库 绿洲 棉花 生长 期 蒸 散发 2009 
年 为 237.86~843.95 mm (均值 为 686.80 mm, 日 均 
蒸 散 发 为 3.71 mm) , 蒸 散发 最 大 值 为 花 铃 期 ,其 次 
为 吐 象 期 .一 期 . 昔 期 和 出 昔 ( 图 4k)。 茸 散发 高 
值 主要 分 布 在 渭 库 绿洲 内 部 河流 周围 及 棉田 小 规 
模 种 植 区域 ( 图 4a~e)。2020 年 棉花 生长 期 蒸 散发 
为 493.13~996.3 mm (均值 为 738.66 mm, 日 均 蒸 散 
发 为 3.99 mm) ,其 中 薰 散发 最 大 值 为 花 铃 期 ,其 次 
为 吐 架 期 . 苗 期 . 蕾 期 .出 苗 ( 图 4)。 蒜 散发 高 值 主 
要 在 绿洲 内 部 及 塔里木 河北 岸 边缘 (图 4f~j) ,由 于 
地 势 低洼 土壤 含水 量 较 高 又 毗邻 沙漠 边缘 , 
偏 高 , 低 值 主要 分 布 在 渭 库 绿洲 东北 和 西南 边缘 ， 
由 于 棉花 集中 连 片 种 植 ,使 用 滴 管 技术 降低 水 分 
蒸发 。 整 体 上 ,12 a 间 棉花 生长 期 蒸 散发 2020 年 
高 于 2009 年 ,其 增长 率 为 7.02% ,尤其 在 花 铃 期 和 
叶 妹 期 蒸 散发 均值 明显 高 于 2009 年 ,但 苗 期 蒸 散 
发 均值 低 于 2009 年 ,然而 取 值 范围 变 大 ,表明 该 
阶段 渭 库 绿洲 气候 差异 性 较为 明显 ,其 余生 长 阶 
段 蒸 散 发 均值 较为 稳定 (图 4k~1)。 随 着 全 球 气 温 
上 升 , 表 明 该 区 域 蒸 散发 也 呈 上 升 趋势 ,棉花 蒸腾 
作用 加 剧 。 
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图 2 2009 一 2020 年 棉花 产量 统计 
Fig.2 Statistics on cotton yield from 2009 to 2020 
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3.3 棉花 水 分 生产 率 时 空 分 析 

3.3.1 和 棉花 水 分 生产 率 时 室 变 化 分 析 基于 像 元 统 
计 水 分 生产 率 频 率 统计 ,可 以 精准 量化 逐 像 元 棉花 
水 分 生产 率 情况 ,由 图 5 可知,2009 年 棉花 水 分 生产 
率 变化 范围 为 0~0.40 kg.m-3( 均 值 为 0.21 kg*m)， 
2020 年 棉花 水 分 生产 率 变 化 范围 为 0~0.45 kgm” 
(均值 为 0.25 kgem™°) ,12 a 间 增长 率 为 16%。2009 
年 棉花 水 分 生产 率 在 0.18~0.25 kg.m-: 之 间 , 有 明显 
的 一 个 峰值 期 ,在 0.21 kgm 对 应 的 像 元 数 最 多 ,达到 
14632 个 。2020 年 棉花 水 分 生产 率 在 0.17~0.33 kg + m” 
之 间 , 有 明显 的 一 个 峰值 期 ,在 0.28 kgm” 对 应 的 像 
元 数 最 多 ,达到 26478 个 ,2009 年 峰值 比较 集中 ， 
2020 年 向 横 轴 右 侧 移动 ,高 值 区 间 增 大 且 较 为 分 
散 。12 a 间 水 分 生产 率 较 高 的 区 域 增多 ,棉花 水 分 
生产 率 整 体 提高 ,布局 更 加 合理 。 将 水 分 生产 率 累 
积 频 率 为 95% 的 像 元 值 统计 作为 临界 值 ,2009 年 对 


(d) 滑 库 绿洲 齐 满 镇 
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2020 年 棉花 产量 /kg'hm- 
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注 : 图 a~e 不 同 颜色 的 点 ,代表 不 同年 份 的 产量 重心 。 
图 3 棉花 产量 标准 椭圆 差 与 重心 迁移 变化 


Fig.3 Standard deviation ellipse and gravity movement of cotton yield 


应 的 值 为 0.37 kgm’, 20204% 0.41 kg.m” ,大 于 该 
临界 值 的 为 水 分 生产 率 较 高 区 域 ,该 区 域 在 2009 年 
主要 集中 分 布 在 渭 干 乡 ,2020 年 分 布 于 玉 奇 喀 特 
乡 ,2 个 乡 均 位 于 新 和 县 且 两 乡 毗 邻 , 土 壤 以 残余 内 
陆 盐 士 为 主 ,灌渠 设施 密集 ,其 中 棉花 呈现 连 片 化 
和 规模 化 种 植 。 

为 准确 分 析 水 分 生产 率 在 空间 的 聚集 程度 和 
演变 趋势 ,以 县 级 行政 和 绿洲 为 单位 ,建立 棉花 种 
植 区 标准 差 椭 圆 和 水 分 生产 率 重 心 点 。2009 年 和 
2020 年 新 和 县 水 分 生产 率 重 心 点 分 布 位 于 渭 干 乡 
和 玉 奇 喀 特 乡 (图 6a) ,标准 差 椭圆 的 长 轴 持 续 增 长 
目 短 轴 持 续 缩短 (图 6e), 自 2009 年 (82.25°E,41.27°N) 
移动 至 2020 年 (82.21°E,41.25°N), 自 东北 向 西南 方 
向 平均 移动 速度 为 888 ma', 由 于 棉花 种 植 面积 的 
扩张 和 产量 的 增加 ,水 分 生产 率 重 心 趋势 逐渐 向 绿 
洲 外 缘 延 伸 。2009 年 和 2020 年 沙 雅 县 水 分 生产 率 
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Fig. 4 Evapotranspiration during cotton growing 


重心 点 均 位 于 古 勒 巴 格 乡 (图 6b), 自 2009 年 (82.53°E， 
41.11°N) 移 动 至 2020 年 (82.52°E,41.12°N) ,重心 点 自 
东南 向 西北 平均 移动 速度 为 320 m-a', 12 a 间 变 
化 趋势 较 小 ,由 于 靠近 塔里木 河岸 和 沙 雅 县 南部 沙 
漠 边 缘 棉花 生长 期 耗 水 量 较 大 ,使 得 水 分 生产 率 重 
心 向 绿洲 内 部 移动 。2009 年 和 2020 年 库 车 市 水 分 
生产 率 重 心 均 位 于 墩 阔 坦 镇 (图 6c ) ,生产 率 标准 差 
椭圆 的 长 轴 缩 短 而 短 轴 增 长 ,其 水 分 生产 率 重心 点 
从 2009 年 (83.05°EE,41.28°N) 移 动 至 2020 年 (83.12°E， 
41.29°N) , 自 西向 东平 均 移 动 速度 为 1071 mar! 
12 a 间 水 分 生产 率 重心 移动 速度 依次 为 : 库 车 市 > 新 
和 县 > 沙 雅 县 ,表明 库 车 市 棉花 水 分 生产 率 以 最 快 
速度 自 西向 绿洲 东部 边缘 移动 ,棉花 产量 和 蒸 散发 
的 变化 直接 影响 水 分 生产 率 重 心 的 迁移 ,同时 说 明 
库 车 市 棉花 种 植 面积 和 产量 重心 癌 绿 洲 边缘 延伸 
面积 最 大 。 

总 体 上 , 渭 库 绿洲 水 分 生产 率 重 心 由 2009 年 


(82.51°E,41.22°N) 移 动 至 2020 年 (82.52°E,41.21°N)， 
自 东 北向 西南 移动 距离 为 1832 m( 图 6d), 年 均 移动 
速度 为 166.55m.a ,表明 绿洲 棉花 水 分 生产 率 呈 现 
扩张 趋势 且 东 西方 向 大 于 南北 方向 。 渭 库 绿洲 水 
分 生产 率 标准 差 椭圆 的 长 轴 持 续 增 长 同时 短 轴 缩 
短 (图 6f) ,长 轴 短 轴 的 共同 作用 下 离心 率 逐 渐 增 长 ， 
同时 标准 差 椭 圆 的 面积 由 2009 年 的 3487.72 kim? X 
渐 增 加 到 2020 年 的 4122.78 km ,县 级 尺度 (图 6e) 和 
流域 尺度 (图 6f) 标 准 差 椭圆 面积 均 增加 ,表明 12 a 
间 渭 库 绿 洲 柳 花 水 分 生产 率 空 间 分 布 方向 趋势 增 
强 ,空间 格局 趋向 集聚 化 。 

3.3.2 棉花 水 分 生产 率 局 部 宝 间 自 相关 分 析 以 $00 mx 
500 m 格 网 数据 为 基础 ,对 棉花 水 分 生产 率 进行 LI- 
SA ZR E YE Fil Moran’s 1 指数 分 析 , 其 中 2009 年 和 
2020 年 棉花 水 分 生产 率 的 Moran’s 1 指数 分 别 为 
0.349 和 0.496, 表明 水 分 生产 率 具有 显著 的 空间 自 
相关 。LISA 聚 类 图 是 对 Moran’s 1 指数 中 通过 了 显 
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1.4 | 一 2020 年 À 主要 位 于 新 和 县 西 侧 绿 洲 边 缘 和 沙 雅 县 古 勒 巴 乡 ， 
al 其 余 高 值 聚 集 区 呈 零 散 分 布 ,该 区 域 棉花 种 植 较为 
ool 密集 , 且 产 量 较 高 ,同时 耗 水 量 低 ,使 得 水 分 生产 率 
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图 5 2009 年 与 2020 年 水 分 生产 率 频率 分 布 
Fig. 5 Frequency distribution of water productivity 
in 2009 and 2020 


著 性 检验 的 区 域 单元 的 地 理 表达 , 渭 库 绿洲 正 相关 
聚集 (高 -高 聚集 或 低 - 低 聚集 ) 较 为 明显 ,聚集 程度 
呈现 上 升 趋势 ,异常 区 (高 - 低 聚 集 或 低 - 高 聚集 ) 呈 


(a) 新 和 县 


N | |O DHE 


玉 奇 喀 特 乡 


高 而 聚集 (图 7a) 。 低 值 聚 集 区 零散 条 带 状 分 布 于 
绿洲 中 部 且 在 沙 雅 县 境内 最 多 ,说 明 棉 花 产量 较 低 
且 耗 水 量 大 ,导致 水 分 生产 率 低 。 异 常 区 零散 均匀 
分 布 于 整个 绿洲 。2020 年 正 相 关 程 度 逐 渐 上 升 且 向 
绿洲 边缘 扩展 ,空间 分 布 基本 上 与 2009 年 一 致 ,但 低 
值 在 绿洲 边缘 和 塔里木 河北 岸 增加 ,表明 该 地 水 分 
生产 率 较 低 ,由 于 新 开垦 的 土地 土壤 贫 盗 ,没有 形成 
大 规模 化 种 植 ,产量 较 低 (图 7b)。 以 上 分 析 得 到 ,高 
值 区 主要 分 布 于 沙 雅 县 新 盟 农 场 和 新 和 县 桑 塔 木 农 
场 , 所 以 在 一 定 程度 上 说 明 ,农场 规模 化 种 植 增产 和 
科学 化 灌溉 ,促使 棉花 产量 和 水 分 生产 率 提高 。 
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图 6 棉花 水 分 生产 率 空 间 分 布 及 重心 迁移 


Fig.6 Spatial distribution and gravity movement of cotton water productivity 
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图 7 水 分 生产 率 LISA 聚集 图 
Fig.7 LISA aggregation diagram of water productivity 


立 人 研究 区 遥感 模型 ,分 析 棉花 水 分 生产 率 的 时 空 分 
布 ,为 渭 库 绿洲 提高 农业 用 水 效率 和 缓解 水 资源 供 
需 矛 盾 提 供 新 思路 。 同 时 ,为 今后 新 疆 小 麦 .玉米 
等 主要 作物 水 分 生产 率 遥 感 估算 提供 基础 。 

在 棉花 产量 遥感 估算 中 ,本 文 使 用 NDVI 指 数 
法 和 统计 数据 空间 化 估算 棉花 产量 ,该 方法 比较 通 
用 且 简 单 ,得 到 的 结果 与 统计 产量 数据 一 致 ,能 
满足 像 元 级 产量 推算 日 效果 较 好 。 以 往 较 多 研究 六 
使 用 增强 植被 指数 (EVT) 建 立 估 产 模 型 ,今后 应 与 
该 指数 结合 对 棉花 面积 提取 和 佑 产 ,进一步 提高 估 
产 精 度 。 从 渭 库 绿洲 棉花 生长 期 蒸 散发 来 看 ,本 文 
使 用 的 SEBAL 模型 遥感 估算 棍 花 生长 期 蒸 散发 较 
好 的 结果 ,在 精度 分 析 中 ,除了 使 用 气象 站 点 数据 
外 ,还 使 用 了 MOD16A2 数据 进一步 验证 ,为 估算 结 
果 提 供 有 利 的 数据 支持 和 可 靠 性 。 针 对 干旱 区 作 
物 蒸 散 发 遥感 估算 的 模型 较 多 ,如 周 柳 萍 “使 用 
SSEB 模型 Zhou 等 “使 用 SEBS 模 型 于 兵 呈 使 用 
HTEM-ABL 模 型 ,今后 研究 中 应 结合 多 种 模型 对 作 
物 生 长 期 的 蒸 散发 估算 。 充 分 定量 化 利用 空间 自 
相关 ,标准 差 椭圆 与 重心 转移 等 空间 分 析 方法 进行 
棉花 水 分 生产 率 时空 分 布 。 本 研究 遥感 数据 估算 
了 棉花 水 分 生产 率 起 始 和 最 终 2 a, 今 后 可 以 补充 并 
加 强 连 续 性 数据 ,如 使 用 高 空间 分 辨 率 Landsat 数据 
和 高 时 间 分 辩 率 MODIS 影 像 时 空 融 合 的 方法 "提高 
水 分 生产 率 时 间 序 列 的 分 析 。 此 外 ,以 往 研 究 中 缺 
乏 新 疆 棉 花 水 分 生产 率 的 研究 ,本 研究 拓展 并 利用 
GEE 平 台 优 势 和 遥感 数据 ,大 面积 对 棉花 水 分 生产 


率 进 行 县 域 尺度 和 流域 尺度 时 空 分 析 与 研究 。 影 
响 水 分 生产 率 的 主要 因素 为 产量 和 蒸 散发 ,12 a 间 
产量 的 增长 速度 超过 了 蒸 散发 的 上 升 速度 ,使 得 水 
分 生产 率 提高 。 

5 结论 

(1) 2008 一 2020 年 棉花 种 植 面 积 整体 呈 增 加 趋 
势 ,棉花 种 植 面积 逐年 向 绿洲 边缘 扩展 。 同 时 , 棉 
花 产 量 也 呈 增 加 趋势 ,从 2009 年 的 1610.10 kg:hm”? 
增长 到 2020 年 的 1855.05 kg:hm ,产量 增长 率 为 
13.20% ,平均 每 年 增长 1.1%, 新 和 县 的 棉花 产量 最 
高 , 库 车 市 产量 最 低 。 产 量 重心 整体 自 西向 东 移 动 
2485 mo 

(2) 渭 库 绿洲 棉花 生长 期 2009 年 区 散发 均值 
为 686.80 mm, 日 均值 为 3.71 mm,2020 年 为 738.66 mm, 
日 均值 为 3.99 mm, 整 体 呈 上 升 趋势 ,增长 率 为 
7.02% ,2009 年 和 2020 年 蒸 散 发 最 大 值 主要 分 布 在 
绿洲 内 部 与 塔里木 河北 部 边缘 。 

(3) 2009 年 水 分 生产 率 均 值 为 0.21 kg.m ，， 
2020 年 为 0.25 kg-m?, 12 a 间 水 分 生产 率 均值 增长 
率 为 16%。 在 空间 上 , 渭 库 绿洲 水 分 生产 率 重心 由 
2009 年 (82.51"F ,41.22°N ) #3) Æ 2020 4 (82.52°K, 
41.21°N) ,在 红旗 镇 自 东 北向 西南 移动 了 1832 m ,年 
均 移动 速度 为 152.67 mea! ,棉花 水 分 生产 率 呈 现 扩 
张 趋势 且 东 西方 向 大 于 南北 方向 ,棉花 水 分 生产 率 
空间 分 布 方向 趋势 增强 ,空间 格局 趋向 集聚 化 。 
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Remote sensing estimation of cotton water productivity in Ugan-Kuqa River 
Oasis based on Google Earth Engine 


HE Xugang'™™, Mamat SAWUT'™, SHENG Yanfang'’’, LI Rongpeng 
(1. College of Geography and Remote Sensing Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 2. Xinjiang Key 
Laboratory of Oasis Ecology, Xinjiang University, Urumqi 830046, Xinjiang, China; 3. Key Laboratory of Smart City and 


Environment Modelling of Higher Education Institute, Xinjiang University, Urumqi 830046, Xinjiang, China) 


Abstract: An accurate and quantitative evaluation of crop water productivity is the basis for the improvement of 
crop yield in arid regions and is of great significance to alleviate water shortage and drive sustainable agricultural 
development. We have considered the Ugan-Kuqa River Oasis on the north bank of the Tarim Basin of China as 
the subject area, and the watershed SEBAL remote sensing evapotranspiration model and cotton distribution iden- 
tification model are established employing the Google Earth Engine cloud platform from 2009 to 2020. The pro- 
duction estimation model was used to evaluate the moisture productivity of cotton in the watershed. The follow- 
ing results were observed. (1) The cotton production in the Ugan-Kuqa River Oasis ranges from 1610.10 kg-hm” 
in 2009 to 1855.05 kg: hm in 2020, with a growth rate of 13.20%. The cotton planting area extends to the edge 
of oasis year by year, and the center of gravity of cotton production moves 2485 m from west to east. (2) The av- 
erage evapotranspiration (ET) value during the cotton growth period was 686.80 mm in 2009 and was 738.66 mm 
in 2020, indicating an overall upward trend, with a growth rate of 7.02%. The maximum value of ET during the 
cotton growth period was observed during the flowering and boll-opening stages. The higher ET value is mainly 
distributed in the interior of the oasis and at the edge of the north bank of the Tarim River. (3) The average water 
productivity in 2009 and 2020 was 0.21 kg:m and 0.25 kg-m , respectively, indicating a growth rate of 16% in 
12 years. In terms of space, the gravity center of water productivity in the Ugan-Kuqa River Oasis moved 1832 m 
from the northeast to the southwest of Hongqi Town, with an average annual moving speed of 152.67 m-a '. The 
water productivity of oasis cotton showed an expanding trend in the east-west direction compared with that in the 
north-south direction. The trend of spatial distribution direction had been enhanced, and the spatial pattern tended 
to be agglomerated. (4) After 12 years, the increased yield rate exceeded the increased evapotranspiration rate, 
promoting the increase of water productivity. Moreover, water productivity is closely related to the cotton plant- 
ing area and reasonable water irrigation technology. High water productivity is mainly distributed in Xinken 
Farm in Xayar County and Sangtamu Farm in Xinhe County, respectively. The large-scale planting and inten- 
sive management of farms promoted and improved cotton production, stable distribution, and efficient use of agri- 
cultural water resources. 

Key words: Google Earth Engine (GEE); evapotranspiration; cotton extraction; cotton yield; water productivi- 


ty; Ugan-Kuqa River Oasis 


